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Stereoscopic Observations of the Blazar 3C 66A with the MAGIC Telescopes 
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Abstract: 3C 66A is an intermediate-frequency peaked BL Lacertae object of uncertain redshift. We report recent 
observations of the region around the blazar with the MAGIC telescopes. The source was observed and detected in 
2009 December and 2010 lanuary, in 2.3 h of good quality data. The signal could clearly be assigned to the blazar 3C 
66A, statistically and systematically rejecting the nearby radio galaxy 3C 66B as a possible origin of the gamma-ray 
signal by 3.6 standard deviations. The derived integral flux above 100 GeV is 8.3 % of the Crab Nebula flux, and the 
energy spectrum is reproduced by a power law of photon index 3.64 ± 0.39(stat.) ± 0.25(sys.). Within the errors, 
this is compatible with the spectrum derived by VERITAS in 2009. From the spectra corrected for absorption by the 
extragalactic background light, we only find small differences between the four modellings that we applied, and constrain 
the redshift of the blazar to z < 0.68. 
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1 Introduction 

3C 66A is a very high energy (VHE) gamma-ray blazar 
classified as a BL Lac object lE6l . or an intermediate- 
frequency peaked BL Lac object (IBL, 11321 ). It has an 
uncertain redshift reported to be 0.444 |29] EJ) or 0.321 
[40], estimated from a single spectral line, or a marginally 
resolved host galaxy respectively. Besides that, the exist- 
ing upper and lower limits to the redshift are > 0.096 [15], 
< 0.44 (331 and < 0.58 ATI . 

In the VHE band, the source was first claimed by the 
Crimean Astrophysical Observatory above 900 GeV with 
an integral flux of (3 ± 1) x 10" 11 cm^s -1 E3- Later 
observations by HEGRA and Whipple reported upper lim- 
its Bl EH. while STACEE found a hint of a signal at 
a significance level of 2.2 fl2l . More recent VERITAS 
observations of 3C 66A took place in 2007 and 2008, 
for a total of 32.8 hours, and resulted in a clear detec- 
tion in VHE gamma rays 0. The derived energy spec- 
trum was compatible with a power law of photon index 
r = 4.1 ± 0.4 sta t ± 0.6 sys and an integral flux above 
200 GeV of(1.3±0.1)xl0- 11 cm- 2 s- 1 (6% of the Crab 
Nebula flux). 

In the GeV band, the gamma-ray emission spot 
3EGJ0222+4253 measured by EGRET was associated to 
3C 66A, although an influence by the nearby pulsar PSR 
J02 18+4232 could not be excluded EUHH. Fermi/LAT 
has monitored 3C 66A since 2008 August, covering the 
latter part of the VERITAS observation. According to Q], 



which reported the first 5.5 months of Fermi/LAT data of 
3C 66A, it showed a significant flux variability (a factor 
of 5-6 between the highest and lowest fluxes). In the new 
multiwavelength study published in [2], photon indices for 
the dark period (1.9 ± 0.1 s t a t ± 0.1 sys ) and flare period 
(1.8 ± 0.1 s tat ± 0.1 sys ) were estimated. In combination 
with the VERITAS spectrum, this indicates a softening of 
the spectrum above 100 GeV. 

From MAGIC observations of that sky region in 2007, 
we reported a significant VHE gamma-ray signal cen- 
tered at 2 h 23 m 12 s , 43°0'7". This excess (named 
MAGIC J0223+430) coincides within uncertainties with 
the position of a nearby, Fanaroff-Riley-I (FRI) type galaxy 
3C 66B (361 z = 0.0215,]. Still, judging from the sky- 
plot alone, and taking into account statistical and sys- 
tematic errors, the probability of the emission to origi- 
nate from 3C 66A was 14.6 %. The energy spectrum of 
MAGIC J0223+430 was compatible with a power-law with 
an index of T = 3.1±0.3. The integral flux above 150 GeV 
corresponded to (7.3±1.5) x 10~ 12 cm" 2 s" 1 (2.2 % of the 
Crab Nebula flux). According to [37], the radio galaxy is 
also a plausible source of VHE gamma-ray radiation. Also, 
the recent MAGIC detection of IC 3 1 [ 8 ] , a radio galaxy at 
a very similar redshift (z — 0.0189) indicates that 3C 66B 
is in principle a feasible object to explain all or part of the 
MAGIC detection from 2007. 
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2 Data set and analysis methods 

In August 2009, 3C 66A went into an optical high state 
which was reported by the Tuorla blazar monitoring pro- 
era This outburst triggered new observations by the 
MAGIC telescopes, located on the Canary Island of La 
Palma (28.8° N, 17.8° W, 2220 m a.s.L). The two 17 m 
diameter telescopes use the atmospheric Cherenkov imag- 
ing technique and allow for gamma-ray measurements at a 
threshold as low as 50 GeV in normal trigger mode. 

We observed 3C 66A in several time slots between 2009 
September and 2010 January. However, the "starguider" 
CCD cameras that are used to cross-check the telescope 
pointing only became fully applicable to stereo observa- 
tions in early December. To allow for a high-confidence 
statement on the directional origin of the gamma rays, we 
only used data taken after these upgrades, which were 5.6 
hours in total. Furthermore, we had to discard data with 
low event rates, affected by the exceptionally bad weather 
conditions in that winter. Finally, we had 2.3 hours of good 
quality data left after all quality cuts. They were taken on 
6 days between 2009 December 5 and 2010 January 18, 
partly under low intensity moon light conditions. 

The data were taken using the so-called "wobble" method 
ifTBI . in which the pointing direction alternates every 20 
minutes between two positions, offset by ±0.4° in RA from 
the source. These wobble positions were chosen to be cen- 
tered on 3C 66A, but the small distance to 3C 66B (0.01°) 
allows equal judgment for both sources. The data were 
taken at a zenith distance between 13° and 35°. 

The analysis we present here involved only events that were 
triggered by both MAGIC telescopes. The analysis was 
done with the MARS analysis framework l30l . taking ad- 
vantage both of the advanced single-telescope algorithms 
(e.g. [ 1 1 ) and newly developed stereoscopic analysis rou- 
tines. A paper about these stereoscopic analysis methods 
of MAGIC is in preparation. 

In the skymapping procedure of MAGIC, which is crucial 
for the directional statements, first an exposure model is 
calculated in coordinates relative to the viewing direction. 
From that, a background expectation distribution is sam- 
pled randomly to celestial coordinates. After applying a 
folding with a Gaussian kernel, the background and real 
event distribution are compared using equation (17) of ||251 
as a test statistic. The resulting null-hypothesis distribution 
of this test statistic is mostly resembling a gaussian distri- 
bution, but in some cases has poissonian components as 
well, therefore the detection significance is always taken 
from the (unsmeared, and unmodeled) (^-distributions. 
The performance of the analysis software was checked with 
contemporaneous Crab Nebula data and simulations. The 
achieved angular resolution^ was 0.1° at 100 GeV and ap- 
proaching 0.065° at higher energies. 

The systematic uncertainty on the direction reconstruction 
is a product of the telescope pointing uncertainty and pos- 
sible biases that occur in the reconstruction algorithms in 
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Figure 1 : MAGIC significance test statistic skymap of the 
region around 3C 66A/B for events with reconstructed en- 
ergies above 100 GeV. 

the presence of e.g. star light inhomogeneities or electronic 
imperfections in the hardware. Both the total pointing devi- 
ation and the telescope pointing precision of MAGIC were 
always monitored over the years lfl3l |7), and along with 
studies of contemporary stereo data of known direction 
lead to an estimate of the maximal systematic stereoscopic 
pointing uncertainty of 0.025°. 

The publicly accessible Fermi/LAT dat£0 were ana- 
lyzed using the public software package LAT Science 
Tools v9.15.2, including the Instrument Response File 
P6_V3_DIFFUSE, and galactic, extragalactic and instru- 
mental background models. 

3 Results 

A a skymap of the observed region above 100 GeV is 
shown in Fig.Q] To calculate a conservative detection sig- 
nificance, we investigated the distribution of squared an- 
gular distances (8 2 ) between photon directions and the as- 
sumed source position. The expected background is ex- 
tracted from the opposite side of the field of view, at the 
same distance from the pointing direction. Comparing the 
events at the source position with this expectation we find 
a significance of 5.2 a (see Fig. |2j. 

Unlike in the previous 2007 detection of MAGIC, the emis- 
sion peak this time is clearly located on top of 3C 66A. The 
fitted center of gravity of the excess (small black square in 
Figure[Qi is at a distance of 0.010° ± 0.023° tat ± 0.025° ys 
from 3C 66A, and 0.108° ± 0.023° tat ± 0.025° ys from 
3C 66B. While being compatible with the blazar 3C 66A, 
the statistical rejection power for the emission to emerge 
from the radio galaxy 3C 66B corresponds to 4.6 standard 
deviations. Even considering the unlikely case of a system- 

1 . http://users.utu.fi/kani/lm/index.html 

2. defined as the a of a two-dimensional Gaussian function 

3. http://fermi.gsfc.nasa.gov/ 
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Figure 2: 8 2 distribution with respect to the position of 
3C 66A (8 2 ) for events with reconstructed energies above 
100 GcV. The OFF data are taken from three positions that 
are symmetrical with respect to the telescope pointing di- 
rections. 

atic offset exactly towards the blazar, the rejection signifi- 
cance of 3C 66B is 3.6 a. 

To derive an energy spectrum, we compared four differ- 
ent unfolding algorithms [5] which correct for efficiency, 
smearing and biasing effects in the energy response of the 
detector. Among these, also the so-called forward unfold- 
ing was tested, which essentially is a fit with correlations 
defined by the response matrix. With all unfolding meth- 
ods, we found that the data are well compatible with a 
power law of the form 



— = K ( - 

dE 200 I 200 GeV 



(1) 



Here the photon index is T — 3.64 ± 0.39 s t a t ± 0.25 sys 
and the flux constant 200 GeV of A' 2 oo = 9.6 ± 2.5 sta t ± 
3.4 sys x 10 -11 cm" 2 s" 1 TeV" 1 . The integral flux above 
100 GeV is equivalent to (4.5 ± 1.1) x lO" 11 cm" 2 s" 1 
(8.3 % Crab Nebula flux). The parameters and statistical 
errors are taken from the forward unfolding, while the sys- 
tematic errors reflect the variations among the other unfold- 
ing algorithms, plus several additional uncertainties dis- 
cussed in |6|. The systematic flux uncertainties add up to 
about 36 % in total. Figure [3] displays the spectral points 
derived using the Tikhonov unfolding method [38], and the 
function we fitted through forward unfolding. 

We also analyzed the Fermi data from the same time period. 
We found the flux variability to be insignificant in a week- 
to-week lightcurve. Given the statistical uncertainties of 
the lightcurve, we would be sensitive on 3 a level to flux 
variations of 60 % or greater, and conclude the variability in 
the days we observed must be less than that. The averaged 
flux above 200 MeV is comparable to the one found in [TJ, 
and less than the flux level seen in 2008 October. A single 
power-law model can reproduce the source spectrum, and 
the photon index is again compatible with the one found in 
Q, indicating no significant change in the overall spectral 
shape. 
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Figure 3: Differential energy spectra of 3C 66 A in the pe- 
riod of 2009 Dec and 2010 Jan. The light-shaded area in- 
dicates the 1 <7 range of the observed unfolded power law 
spectrum gained by forward unfolding, the crosses are from 
the unfolding after [38] for comparison. The dark-shaded 
area is the spread of the EBL-corrected, mean flux values 
obtained by the four applied EBL models, assuming the 
redshift of z=0.444. The VERITAS (observed) spectrum 
after [3| is shown for comparison. 

4 Discussion and Conclusions 

MAGIC detected a VHE gamma-ray signal from 3C 66A in 
the period between December 2009 and January 2010, and 
during an optical active state of 3C 66 A J9J. We rule out 
the emission to come from 3C 66B at a confidence level 
of 3.6 a. In our previous detection from 2007 we could 
not significantly reject any of the two sources, so the new 
detection does not contradict the previous conclusion that 
the signal in 2007 might have been emerging from the radio 
galaxy. In particular, both objects are not only viable VHE 
emitters, but both would be likely to be variable sources. 

Given the strong flux we detected in 3C 66A, we further- 
more conclude that 3C 66A might have to be in a low flux 
state in order not to outshine the comparably weak emis- 
sion from 3C 66B, assuming it is responsible for the flux 
level MAGIC detected in 2007. 

The energy spectrum we unfolded is softer than in the pre- 
vious MAGIC detection (T = 3.10 ± 0.31 sta t ± 0.2 sys ), 
and at the same time compatible with the VERITAS spec- 
trum of 3C 66A. Compared to VERITAS, MAGIC has a 
lower threshold and the spectrum is extending to well be- 
low 100 GeV. The flux level of 8.3 % Crab Nebula flux 
is similar to the one reported by VERITAS (6 %), and sig- 
nificantly higher than in the previous MAGIC observation 
(2.2 %). 

The VHE gamma rays produced at the source can be ab- 
sorbed in the intergalactic space by pair production with 
the UV to infrared photons of extragalactic background 
light (EBL) Il34l [191 . The amount of absorption depends 
on the energy and redshift, and can be corrected for in 
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the data, assuming a given modeling of the EBL density. 
This kind of de-absorbed spectrum can be regarded as the 
source-intrinsic spectrum that we would measure if there 
were no EBL. To derive such a de-absorbed spectrum, we 
tested several state of the art EBL models, namely ifTTl . 
the fiducial model in (H, E2 and OH. The EBL correc- 
tions were applied in the unfolding procedure, because it 
needs the covariance matrix to correctly calculate the er- 
rors. The spread of the differential, de-absorbed flux spec- 
tra, obtained with the four models and assuming the red- 
shift of z = 0.444, is shown as the dark shaded area in 
Figure[3] The de-absorbed photon indices for the four EBL 
modelings we used are 2.57 ± 0.68 Q7), 2.61 ± 0.67 d, 
2.59 ± 0.68 CH and 2.37 ± 0.70 (22|. We conclude that 
the differences between the de-absorbed spectra are very 
small, reflecting the fact that also the predicted EBL shapes 
and densities are very similar. 

Following the predicions of most VHE emission models, 
the de-absorbed spectrum of a blazar is not expected to be 
concave, i.e. rising towards higher energies. Two ways of 
testing this are to comparing the various points of our own 
spectrum, or to compare the points with the Fermi photon 
index. The fact that we find our spectrum neither signifi- 
cantly concave, nor harder than in Fermi, suggests that the 
assumed redshift of z = 0.444 does not contradict our ob- 
servations. We investigated the plausibility of the redshift 
quantitatively, assuming that the intrinsic spectrum is not 
expected to be exponentially rising, and thus have a pile- 
up, at highest energies. This common method was previ- 
ously used and described for example in l27ll28l . Using the 
IfTTl model and the likelihood ratio test between the "power 
law" and "power law + pile-up" hypotheses, as described 
in the reference, we derive an upper limit on the redshift of 
z < 0.68. 

Achieving these results in only 2.3 h of observations 
demonstrates the striking advantages of the MAGIC stereo- 
scopic system with respect to its monoscopic era. Further 
MAGIC and other gamma-ray observations of this region 
can provide interesting information about the IBL type BL 
Lac object 3C 66A, and, during low flux periods of that, 
also the FRI type galaxy 3C 66B. 
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